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Abstract 

Accretion disks with an additional mass input on the disk surface 
from environment are considered in the limit of low mass input rate, 
i.e. when the accretion flow remains keplerian. Due to dissipation 
of kinetic energy of the infalling gas disk temperature increases, and 
can deviate significantly from standard temperature of reprocessing 
accretion disks in their outer regions. This increase in temperature 
produces an excess of emission in a long-wavelength range of the disk 
spectrum. An illustrative example of the spectrum of a weakly mass- 
loaded disk (the surface mass input rate Mi oa d — 1-5 X 10~ 8 M© yr _1 ) 
reasonably reproducing FIR excess observed in CS Chameleon is given. 
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1 Introduction 



Dynamics of the multi-phase interstellar medium is dominated by mass 
and momentum exchange between different phases. Ejection of mass into 
diffuse component from a condensed phase changes velocity of the bulk flow, 
increases pressure in diffuse phase, and heats diffuse gas due to conversion 
of kinetic energy of evaporated/destroied clouds into thermal energy of bulk 
medium (Cowie, McKee & Ostriker 1981, Hartquist et al. 1986). As a re- 
sult, loaded flows (i.e. the flows dominated by mass input from a condensed 
phase) can differ qualitatively from the ones in a homogeneous (monophase) 
medium. In a simplest approach mass loaded flows are described by dynami- 
cal equations with source terms. For accretion disks mass loading corresponds 
physically to a continuing mass infall on the disk surface from the residual 
gas envelope. 

Extremely inhomogeneous structure of the diffuse interstellar gas and 
molecular clouds shows self-similarity on the scales from 0.02 (the resolution 
threshold) to 100 pc (see, Falgarone, Phillips & Walker 1991). It seems obvi- 
ous that such self-similarity can extend to much smaller scales corresponding 
to protoplanetary condensations in accretion disks. [Note that the resolu- 
tion limit, 0.02 pc, corresponds to planetary masses of clumps, ~ 10~ 4 M ]. 
One can expect that in such conditions clumpy gas from debries of the par- 
ent molecular cloud will provide a continuing matter infall on the accretion 
protoplanetary disk being already formed. In theoretical treatment of pro- 
toplanetary disks this circumstance is described by the source terms in mass 
and momentum equations (see, e.g. Safronov & Vityazev 1983, Papaloizou & 
Lin 1995). Multiple CI and CO velocity components (Roberge et al. 2000), 
Call absorptions (Barnes, Tobin & Pollard 2000) and emission in Fell, SI, 
Sill, Nill and CI lines (Lecavelier des Etangs et al. 2000), observed in (3 
Pictoris, are seemingly connected with the infall of clumpy gas on the disk. 
Similar detections of the infalling circumstellar gas have been reported for 
Herbig Ae/Be stars (Grinin, Natta & Tambovtseva 1996). One can assume 
that this phenomenon represents late stages of matter infall from a highly 
inhomogeneous primordial circumdisk cloud onto the accretion disk. If this 
is the case, and protoplanetary disks are loaded by infalling and passing 
through them dense clumps, both vertical and radial structure of such disks, 
the accretion rate and the temperature distribution can differ considerably 
from those in homegeneous (monophase) disks. In this paper we analyze 
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properties of steady disks, concentrating mainly on radial temperature dis- 
tribution, in a simple model of weak mass loading, i.e. when mass loading 
does not violate the keplerian flow of the bulk motion. In Section 2.1 we 
describe the model, in Section 2.2 we formulate basic equations describing 
radial structure of thin accretion disks for weak mass loading limit, in Section 
2.3 we describe the solution of energy equation and find radial temperature 
distribution, and discuss possible observational manifestations, the effects 
of angular momentum input associated with mass loading on temperature 
distribution are discussed in Sec. 2.4, Section 3 summarizes the results. 

2 Model and governing equations 

2.1 Model of weakly loaded disks 

As a model of a mass-loaded disk which admits simple analytical de- 
scription we will accept the following: an accretion disk is surrounded by 
a spherical "subsystem" of dense condensations formed of the debries of a 
parent cloud and orbiting around the central star (Fig. 1). This subsystem 
can be in a quasistationary state either because the condensations are nearly 
dissipationless, or an infinite (in practice, very massive) reservoir of gas re- 
plenishes dissipation of mass and energy. In the first case, the characteristic 
life time of the quasistationary orbits determined by interaction of the con- 
densations with the diffuse component, with other condensations and with 
the disk itself, must be much larger than characteristic rotation period of the 
disk. For this to fulfill the following relations must hold, respectively 

iV,<iV c , (iV c )<iV c , N d ^N c , (1) 

where iVj, N c , and are the column densities of the diffuse intercloud 
gas over the whole circumdisk cloud, of a single clump, and of the disk, 
respectively, (N c ) = fN c , f = M C R 2 JR 2 is the covering factor of clumps, R c 
is the clump radius, R, the radius of the circumdisk cloud. These inequalities 
definitely fulfill for planetesimals or comet-like bodies. For gaseous clumps it 
is not so obvious, and one has to assume a large enough density contrast in 
order that the clumps were long-living and orbiting quasistationary. In the 
latter case when the clumps are not sufficiently dense, so that for instance 
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N c ~ Nd and they are absorbed by the disk in a single encounter, the overall 
picture can be kept quasistationary due to replenishment of the clumps from 
the surrounding molecular cloud as shown schematically on Fig. 1. A gas 
reservoir with radius R ~ 3 — 10 pc can support continuous mass infall during 
the time t ~ 1 — 3 Myr with the rate 10~ 7 — 10~ 6 M yr^ 1 if the mean density 
is only n c f v ~ 0.03 — 0.3 cm -3 , where n c is gas density in clumps, f v , the 
volume filling factor of the clumps; optical depth corresponding to the whole 
ensemble of such clumps is as small as r c ~ lO~ 21 n c f v R ~ 0.03. As we will 
see below, mass loading even with one order of magnitude lower mass rate 
can produce observable effects. We will assume that the quasistationarity 
takes place in either way. With these assumptions the subsystem of clumps 
can be treated as stationary. When the clumps intersect the disk they lose 
a fraction of their mass due to stripping under the action of Rayleigh- Taylor 
instability (see detailed discussion in Hartquist et al. 1986, Klein, McKee 
& Colella 1994). The subsystem of clumps is assumed in rest as a whole, 
i.e. (v c (r)) = 0, where the averaging is over the ensemble of clumps at given 
r, such that the angular momentum input into the disk from destroying 
clumps is zero. More general case with v c = v c (r) will be discussed in Sec. 
2.4. 

2.2 Equations of motions 

In this framework the continuity eqation for thin disk is described by the 
equation 

where £ = f*™ pdz, ur = v r dz, (the interrelation T,ur = f*™ pv r dz 
is assumed as well), * = J*™ qdz, q [g cm -3 s" 1 ] is the mass ejection rate 
from destroying clumps per unit volume. [In general, \1/ may represent not 
only mass loading due to destroying cloudletts, but also mass infall onto the 
disk from the envelope (Safronov & Vityazev 1983), the effects of the mass 
transfer stream in cataclismic variable stars, and the mass loss through a 
wind (Papaloizou & Lin 1995)]. 

Naviet-Stokes equation for azimuthal component of a stationary disk with 
zero angular momentum input is written in a standard form (see, Shakura & 
Sunyaev 1973, Pringle 1981) 
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± { ^u R ) = A („ Efl ^) , (3) 

here = is explicitly taken into account, Q = Q(R) is the local angular 
velocity, v is the kinematical viscosity. The change of the angular velocity Q 
in one rotational period Tr = 2n/Q due to mass input is \AQ\ ~ ^/E. From 
requirement that |Af2|/f2 <C 1 one obtains restriction on the mass loading 
\l/ <C Y*Q/2tt, so that the ratio of the mass input rate M p associated with the 
infalling clumps to the standard mass accretion rate M ~ 2-rvYj is 

M„ 1 R 2 Q 1 us 

-J- < ~ -TT^l > L 4 

M 4tt i/ 4vr 
Obviously, the restrictions on the muss input \I/ are rather weak. Formally, 
the change in Q is larger on a much longer accretion time, Tp = R/\ur\ ~ 
R 2 jv (here ur ~ z//i? is explicitly assumed, see Pringle 1981) 

\AQ\~n Mp, (5) 
1 1 M ' 

where Qq is characteristic angular velocity in the absence of mass loading. 
For the radial velocity component we have 

du R ul 1 dp GM ^u R 
UR l[R-R + -pdR + lV = -— (6) 
Assuming that for a weak loading the radial component ur is of the order 
v/R (see, Pringle 1981) we arrive at 



GM r 



n [1 + 0(M- 2 ) +P(R)\, (7) 

where f3 = ^urR 2 /GMT,, M. = us/c s is the local azimuthal Mach number, 
c s , the local sound speed. Equations ©-© differ from the standard equa- 
tions of thin accretion disk (Shakura & Sunyaev 1973, Lyndel-Bell & Pringle 
1974, see also reviews by Pringle 1981, and Papaloizou & Lin 1995) by the 
sources due to mass loading in the r.h.s. of continuity equation (fj) and in 
equation for radial velocity (|6|). The parameter (3 can be estimated as 

aT R M p 

7T 1 D Mq 
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with the assumption that M p /Mq can have low or moderately high values, 
here a is the viscosity parameter (Shakura & Sunyaev 1973), and thus its 
contribution to the r.h.s. of (|7|) is negligible, and in this sense the loading 
will be treated as weak. It is readily seen that approximately (3 ~ 27r\l//fi£, 
and (ID is equivalent to the condition that on one rotational period |Af2| <C 
Q. This means that the infalling settles onto the generic keplerian motion 
quickly, and as long as the mass input from clumpuscules is small enough on 
one rotational period, i.e. (Q) fulfills, it can be treated as a weak loading. 

Defining the accretion rate as M = —27rKHuii (see Pringle 1981) one can 
obtain the following solutions of equations (0) and (||) for a weakly loaded 
(i.e. p < 1) disk 

R 

M = M* - 2tt J ^RdR, (9) 

R 



3tt 



R \ 1/2 



^j^RdR, (10) 

R* 



where M* = — 2it(RYiUr)* is the total mass accretion rate (i.e. the rate on 
the stellar surface); note, that the mass rate provided by the accretion disk 
itself in the absence of mass loading is equal 

oo 

M = M, - 2tt j ^RdR. (11) 
The quantity z/S determines dissipation at a rate 



D(R) = \uH(R^') 2 = ^^[i - (RJR) 1/2 \ - ^ / *(R)MR, (12) 

R* 

and the corresponding luminosity integrated over the disk 



°r 1 CM ' M °r HR ^ 

U D = 4vr / D(R)RdR = - — - 3ttGM / / ^(R')R'dR'. (13) 
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2.3 Energy equation 

One of the most interesting and important consequence of the mass load- 
ing is connected with possible heating of the disk due to conversion of kinetic 
energy of the falling gas into heat. The energy equation in the presence of 
mass loading (see, Cowie, McKee & Ostriker 1981, White & Long 1991), 
integrated over z, is written as 



2 fi 



r ±(7 
2 W 



P(u 2 ) 
P 
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£(T,p), (14) 



where e is specific (per unit mass) thermal energy, the first term in the r.h.s. 
describes the energy input due to inelastic interaction of gas stripped from 
clumps and involved in keplerian motion of the disk, (u 2 ) is the mean square 
relative velocity of clumps and disk (the averaging is over the ensemble of 
clumps), 1Z is gas constant, £ is the radiative cooling rate, we will assume 
C = aT 4 , o is Stefan-Bolzmann constant, the second term is the dissipative 
energy rate ([12]) . In steady case, disk temperature is determined 



aT 4 



3TZ 



x) p{u 2 ) 



2 p 
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i\2 



(15) 



The mean square velocity is (u 2 ) 



11- 



where u c is the velocity disper- 



sion of clupms encountering the disk, for clumps orbiting around the central 
star or falling from external molecular cloud u c ~ u^, so we will assume 
(u 2 ) ~ 2w^; the contribution from radial disk velocity ur is neglected. The 
corresponding contribution to disk luminosity from the first term in the r.h.s. 
of (0) is 



T c 



K f 
6tt— / Tm 
a J 

R* 



, P(^ 2 ) 
p 
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RdR, 



(16) 



and when compared to the viscous luminosity in the absence of mass loading 
it is 



T c 
p 



R 2 m M n 



Mn 



(17) 



In dimensionless units Eq. (|15D reads as 



7 



X 



= 2 (7 -d^|4 (x - -j=) -^1 J ( 18) 

where x = R/R*, ^ assumed in the form \1/ = ty ip(x), & — T/T*, 5 = 
{c*/c M ) 2 



oc 



J = / V^cte, = 7— T*, c 2 M = — ^. (19) 

// it* 
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The problem contains two dimensionless parameters, 5, and A. It is readily 
seen that 8 < I: for i2* = R Q , M = M & and M ~ 10~ 6 M yr" 1 this ratio is 
5 ~ 10~ 3 . For weak loading, i.e. M* £ M , A — 1 is also small, however even 
in this case mass loading can change significantly temperature distribution 
in the disk and its spectrum. In Fig. 2 we show the temperature distribution 
6(x) for the mass- loading function 



rR — -R* 



Ri 



R — R* 



${R) = exp — i + 1 , (20) 



R 



L 



with the peak ip(R) — 1 at R — R* + Rl, where Rl is the scale of the radial 
distribution of matter infall, Rl = LR*, L = 30 and a set of A. 

It is seen that the most important contribution from additional mass input 
comes from regions x ^ L where temperature can increase by factor of 2 even 
for a moderate mass loading A = 2. At larger distances, x ~ 10 — 20 L, the 
temperature profiles for A > 1 approach the standard profile (A = 1), and 
then go lower because at such distances the heat input from kinetic energy 
of clumps (u 2 oc x" 1 , first term in the r.h.s. of Eq. [18]) is smaller than the 
energy required to heat the infalling (presumably cold) gas (second term in 
the l.h.s. of Eq. [18]). Similar changes are seen for more extended in radius 
mass input, L = 100 (Fig. 3). In this case, the surface mass input rate 
for equal A is ~ 3 times lower than for L = 30, as A — 1 oc ^oL, however 
the increase in total luminosity of loaded disks from the area oc L 2 must 



be larger. For mass input in the form (pTJ) and A — l~latx<L/3 the 
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second term in the r.h.s. of ( Jl8| ) dominates and the first and third terms 
contribute only around 0.01 % to 1%, at x ^ L/3 the situation changes and 
main contribution (up to 97%) is from the first tirm in flTE|). At x > 10L 
the dominant contribution is again from the second term; the third term at 
x > L/3 remains always about 30-50 % of the second. For A — 1 ~ 0.1 the 
second term dominates up to x ~ 2L/3 and then at x > 8L, the last term 
always remains as small as 5-10 % of the second term. 

In order to illustrate how the mass loading can change the spectra of 
accretion disks, we show here the spectral energy distribution (SED) for 
optically thick disks seen face-on 

Rout 

Z u = 2n J B v [T(R)]RdR, (21) 

with T(R) = T*9(x), B U [T] being the Planck function, R out , the outer radius 
of the disk, in the models shown in Figs (2)-(5) R out = 500, contribution 
from regions R > 500 is small. Defining dimensionless frequency y as 

hv , , 

y = (22) 

one can reduce ( pT| ) to the form 

K = ZoS(y), (23) 



where 



Wtt^RI 



(24) 



h 2 c 2 

Fig. 4 shows the spectra from disks with temperature distribution shown 
in Fig. 2. It is seen that even for relatively small loading, A = 2, the 
spectrum changes substantially: in maximum it increases by factor of 5, the 
maximum itself is more sharp showing an excess of energy in this range (for 
R = R & , M = M Q and M = 10~ 6 M & yr -1 the frequency at maximum is 
v ~ 1 — 1.5 x 10 13 Hz, or A ~ 10 — 20 /im). For a more extended disk [L = 100) 
the changes in their spectra are more pronounced as shown in Fig. 5: at A = 2 
in maximum it reaches an order of magnitude increase in comparison with 
mass-unloaded disks, the maximum itself shifts to a lower frequencies, and the 
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energy excess in this frequency range is more obvious. These two examples 
show that additional heat input due to mass loading can be responsible (or 
may contribute) for FIR excess observed from circumstellar disks around 
young stars (see for a review Natta, Meyer, Beckwith 1997). In order to 
illustrate such a possibility we compare in Fig. 6 the spectrum expected 
from an optically thick mass-loaded disk seen face-on with the one observed 
in the Chameleon cloud (Natta et al. 1997). For the modelled spectrum 
we assumed the distribution of infalling material of L = 9.2 AU in radius 
(integration was over R = 3L), its accretion mass rate Mq = 2.2 x 10~ 7 M Q 
yr -1 , and the mass loading only A— 1 = 0.08, the central star is of M = 1 M , 
i£* = Rq. Filled squares show the FIR excess observed from CS Cha and 
which cannot be attributed to a standard reprocessing accretion disk; open 
diamonds in optics show presumably the photospheric emission. [The lack, or 
negligibly low level, of emission in optical and NIR range from the accretion 
disk is considered as an indication of a hole in the inner disk (about 0.3 AU 
in size).] 

For a rough estimate of the effects of mass loading on SEDs from accretion 
disks, one can use the following approximate expression for temperature in 
the region where heating due to dissipation of kinetic energy of falling clumps 
(first term in r.h.s. of Eq. pl|) reaches the maximum (at x m ~ \fL with 

{s/x) m ~ e/L) 

*m~[2(7-l)^~] 1/4 , (25) 

which correctly by order of magnitude describes the asymptotics at distances 
x where accretion heating (last two terms in Eq. |18|) and cooling due to mix- 
ing of cold infalling gas (second term in the l.h.s. in Eq. [18]) are unimportant 
- this region corresponds to the plateaus in 6(x) profiles in Figs. 2 and 3. In 
dimensional units T m ~ 0.8(A — 1) 1//4 L _1/,2 T*. For the above model it gives 
T m ~ 100 K, and the corresponding frequency v m ~ 6 x 10 12 Hz. 

2.4 Effects of angular momentum input 

When infalling gas has non-zero angular momentum equation, so that v c ^ ^ 
equation (|3|) is written in the form (Safronov & Vityazev 1983, Papaloizou 
& Lin 1995) 
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A(£i? 3 ^) = A ^ SjR 3 + y( R ) R 2 Vc ^ R ^ (26) 

and the contribution to energy equation is determined by 



R 

2tt 



J ^{R)R 2 v c ^R)dR. (27) 



MqR?VL 

J 

In an extreme case when the angular velocity of the infalling gas approaches 



the keplerian value v c ^(R) = JGM/R energy equation reads as 



A - 1 1 r , , A - 1 1 



x 3 



J tpxdx + — ^ J ipx z/2 dx. (28) 

i X i 

It is readily seen that for the mass input rate in the form (p0|) the contribution 
from the last term is of the same order of magnitude (about 2/3) as of the 
third term, and thus never exceeds 30% for A ~ 1 and 3-6% for A — 1 <C 1. 



3 Conclusions 

Dissipation of kinetic energy of clumpy gas falling onto accretion disks from 
the debries of a parent protostellar molecular cloud significantly alters the 
radial temperature distribution. It can be important in outer parts of a re- 
processing disk, where temperature decreases sufficiently T ~ -R -3 , while the 
heat input rate associated with mass loading is approximately proportional 
to the square of the azimuthal velocity E oc R~ x . Thus, mass loading can 
naturally produce excess of emission in the long-wavelength range of SEDs 
of the accretion disks, and with an appropriate radial distribution of mass 
loading can reproduce observed FIR excess. For a one-peak radial distribu- 
tion of the infall rate localized at R ~ L, the corresponding temperature in 
the hot ring for i?* = 1 R Q , M* = M Q is 

Tm „ 800 (^) (_U ^ ( 6 f° \ K, (29) 

\mJ vi au; v 10 " 6M ©y r / 



ii 



and the wavelength 



Am . 6 (J-f m- 1 " { *> „ m . (30) 

It is readily seen that the extra luminosity produced by the infalling gas 
is proportional to the mass loading rate, L(M p ) ~ aT^L 2 oc M p , and for 
weak mass loading, M p <C Mq, is small as compared to the total luminosity 
of an unloaded disk. However, it can result in an observable effect if the 
loading mass is deposited at large radii, so that the wavelength in the peak 
of emission (|30"D falls into FIR range where the reprocessing disk does not 
contribute singnificantly. This case is shown in Fig. 6 where for mass loading 
M p ~ 0.1M the FIR luminosity of the infalling gas is approximately ten 
times smaller than the optical luminosity of the reprocessing disk. 

The results presented in this paper have only an illustrative meaning, 
and can be directly applied for explanation of spectral features observed in 
protoplanetary disks only when clear evidences of mass loading will be found. 
Nonetheless, they definitely show that mass loading can have an important 
influence on temperature distribution in accretion disks and their spectra. 

I thank S. A. Lamzin, A. V. Lapinov, A. M. Sobolev, I. I. Zinchenko 
for valuable discussions. This work was supported by RFBR (project No 
99-02-16938) and INTAS (project No 1667). This research has made use of 
NASA's Astrophysics Data System Abstract Service. 
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Note added in manuscript 2002 August 8.- After acceptance of this pa- 
per I have been informed about the papers by Falcke H. & Melia F., 1997, 
Astrophys. J., 479, 740, and Coker R. F., Melia F., & Falcke H., 1999, As- 
trophys. J., 523, 642 which deal with accretion disks around massive black 
holes loaded by wind infall with application to Sgr A*, in the framework of 
hydrodynamical equations similar to that used in this paper. Their emission 
spectrum has an excess in the long-wavelength limit due to additional energy 
input connected with the mass-loading similar to that obtained in our paper. 
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Figure 1: Schematic representation of the model: debries of the seed cloud 
falling down on the accretion protoplanetary disk. 
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Figure 2: Temperature distribution 9{x) for the mass-loading function (|20|). 
5 = 10~ 3 , A = 1 (unloaded), 2, 3, 4 from the bottom to top, L = 30. 
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Figure 5: Same as in Fig. 4 for L = 100. 
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Figure 6: SED in CS Cha: open diamonds from ground-based observations 
(presumably represent the spectrum of the stellar photosphere), filled squares 
from ISO and IRAS - Natta, Meyer, & Beckwith (1997), solid line shows 
the spectrum of a mass-loaded disk with A = 1.08, L = 9.2 AU, M = 
2.2 x 10" 7 M yr" 1 . 
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